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Abstract: With the upcoming run of the Large Hadron Collider (LHC) at much higher 
center of mass energies, the search for Beyond Standard Model (BSM) physics will again 
take center stage. New colored particles predicted in many BSM scenarios are expected to 
be produced with large cross sections thus making them interesting prospects as a doorway 
to hints of new physics. We consider the resonant production of such a colored particle, the 
diquark, a particle having the quantum number of two quarks. The diquark can be either 
a scalar or vector. We focus on the vector diquark which has much larger production cross 
section compared to the scalar ones. In this work we calculate the next-to-leading order 
(NLO) QCD corrections to the on-shell vector diquark production at the LHC produced 
through the fusion of two quarks as well as the NLO corrections to its decay width. We 
present full analytic results for the one-loop NLO calculation and do a numerical study to 
show that the NLO corrections can reduce the scale uncertainties in the cross sections which 
can be appreciable and therefore modify the expected search limits for such particles. We also 
use the dijet result from LHC to obtain current limits on the mass and coupling strengths of 
the vector diquarks. 
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1 Introduction 

After the successful running of the Large Hadron Collider (LHC) at CERN with 7 and 8 
TeV center of mass energies, the data released by the two experiments, ATLAS and CMS 
have not only improved on the limits set by the Tevatron experiments on any new physics 
scenario, but has also started giving some insights into the TeV scale. In addition to the 
observation of a scalar resonance at 125 GeV [1, 2] consistent with that of the Standard 
Model (SM) like Higgs boson, the results are also in very good agreement with predictions 
from the SM, with not much deviation. This means that the LHC data is already pushing the 
energy frontier of any Beyond Standard Model (BSM) physics predictions. However with the 
upgraded run of LHC at center of mass energy of 13 TeV and subsequently 14 TeV, the search 
for new physics is expected to be more robust and as envisaged for the LHC run. As expected 
and observed from the previous LHC runs, the data would be most sensitive to the strongly 
interacting sector through production of new colored states. Since the initial states at hadron 
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colliders such as the LHC are colored particles, the most dominant contributions would be 
through new colored resonances. Such colored particles are predicted in many class of BSM 
theories. Resonant s-channel production at LHC can happen for squarks in R-Parity violating 
supersymmetric theories [3], diquarks in super-string inspired Eq grand unification models [4] 
or models with extended gauge symmetries [5-7], color-octet vectors such as axigluons [8, 9] 
and colorons [10-13], models with color-triplet [14], color-sextet [15-17] or color-octet scalars 
[18, 19]. The absence of any such observation in the existing data put strong limits on such 
particle masses, from pair production of such states, or more strongly from resonant searches 
of new physics exchanged in the s-channel [20]. 

These resonant colored states are most likely to decay to two light jets leading to not only 
the modification of the dijet differential cross section at large invariant mass but also show up 
as a bump in its invariant mass distribution. Such a signal will not go unnoticed and will be 
fairly very distinct at large invariant mass values, as the significantly huge QCD background 
falls rapidly for large dijet invariant mass. Both ATLAS and CMS Collaborations have looked 
at the dijet signal and already put strong constraints on the mass of such resonances [21— 
26]. We should however note that the production of such colored particles will be beset 
with significant contributions from QCD corrections, and therefore it becomes important to 
understand how much the leading order (LO) rates might change once these corrections are 
included. One finds that there have been significant efforts in this direction to study the 
next-to-leading order (NLO) QCD effects on production of some of the new colored particles 
[27-29] arising in BSM at the LHC. Here we are interested in particular with particles of 
the “diquark” type which carry non-zero baryon number and couple to a pair of quarks or 
anti-quarks. The fact that LHC being a proton-proton collider will have valence quarks in 
much abundance compared to the anti-quarks, helps in producing the diquark as a resonance 
through qq fusion. A lot of studies carried out at LO exist in the literature for such diquarks 
and their resonant effects in the dijet signal [30-37], pair production of top quarks [38-41] 
and single top quark production at the LHC [42, 43]. The one-loop NLO correction for scalar 
diquark production was considered in Ref [27]. We focus on the case of vector diquarks which 
are either antitriplets or sextets of SU(3)c- Such particles will also be copiously produced 
as s-channel resonances with much larger cross section compared to the scalar ones. Once 
produced, the vector diquark will decay and would thus contribute to the dijet final state or 
to final states involving the third-generation quarks. 

For our study of estimating the NLO corrections to the on-shell production of a vector 
diquark at the LHC, we follow in part the methodology used in Ref [27] to present our results. 
In Sec. 2 we present the formalism and give the basic interaction Lagrangian relevant for our 
study and in Sec. 3 we discuss the on-shell production cross section of the vector diquark, 
and present our calculations and analytic expressions for the NLO QCD results. In Sec. 4 
we give results for the one-loop corrections to the decay width of the vector diquark. In Sec. 
5 we give our numerical results for the NLO cross sections and its dependence on the choice 
of scale for the production of the vector diquark in different channels at the LHC. We also 
consider its effect on the experimental limits for such particles and finally in Sec. 6 we give 
our conclusions with future outlook. Some relevant formulas are collected as an Appendix. 
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2 Formalism 


We are interested in new colored particles that couple to a pair of quarks directly and carry 
exotic baryon number. With the LHC being a proton-proton machine, the initial states com¬ 
prised of the the valence quarks (it, d) would lead to enhanced flux in the parton distributions 
for the collision between a pair of valence quarks such as uu, dd or ud. Any new particle that 


couples to these pairs would carry a baryon number B = — and will be charged under the 
SM color gauge group SU(3)c- Such states are generally referred to as diquarks. These 
colored diquarks can be either color antitriplets or sextets of SU(3)c- We can describe the 
vector diquarks following Ref [43] according to color representation (3, 6) and electric charge 
(4/3, 2/3, 1/3) as V ^ D , V^ D , , where the subscripts 2 U,U, and V in the fields indicate 
their electric charge |Q| of two up type quarks, one up and one down type quark respectively, 
while Nd(= 3 ( 6 )) is the dimension of the antitriplet (sextet) representation. The relevant 
interactions of the quarks with the different vector diquarks is given by the Lagrangian 


r v 

L/ qqD 



\ 2U 

^/ l _l_ § = + 


- 7 =£=v^v c aallx p r Vp b 

y 1 + Oa/3 


+ A VpV^U c aalfJ P T Vp b 


+ h.c. 


( 2 . 1 ) 


where P T = \{ 1 ± 75 ) with r = L,R representing left and right chirality projection operators 
and superscript /a is the Lorentz four vector index. The K J ab are SU(3)c Clebsch-Gordan 
coefficients with the quark color indices a,b = 1 — 3 , and the diquark color index j = 1 — Njj , 
C denotes charge conjugation, while a, (3 are the fermion generation indices. The color factor 
K J ab is symmetric (antisymmetric) under ab for the 6 (3) representation. A more general 
form of the Lagrangian can be found in Ref [ 35 ]. A factor of l/y/2 in the interaction terms 
involving same quark flavors is introduced to keep the expressions for the production cross 
section as well as the decay width same for both different flavor and same flavor cases. To 
calculate the QCD corrections to the diquark production, we also need to know how the 
vector diquark (V/ 4 ) interacts with the gluons, which is given by the Lagrangian 1 * * , 

£ V gdd = (vn - igs vlT$V jv G A ^ ( 2 . 2 ) 

where, 

vr = D§vr - B^V? ( 2 . 3 ) 

G% = d,G A - d u G A + g s f ABC G*G c u ( 2 . 4 ) 
= Sijdf, - ig s G A gT A j. ( 2 . 5 ) 


The indices i and j again run from 1 —> Nd, where Nd is the dimension of the diquark 
representation. The index A runs from 1 —>• 8 and Tfj are the SU(3)c generators in the 

1 There may exist anomalous terms in the Lagrangian allowed by gauge invariance, similar to that for vector 

leptoquarks [44]. For simplicity, we have neglected such anomalous contributions in the gluon-diquark-diquark 

interaction. 
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diquark representation. Note that we have suppressed the electric charge index (2 U,U,V) 
for the diquark as we are interested only in the QCD corrections. The Feynman rules for 
three-point vertices involving vector diquark are given in Appendix A. 

The diquark can couple to the initial state valence partons coming from both the protons, 
and the production of the diquark would get significant enhancement due to the large flux 
of the valence quarks in the proton. Therefore the production rates are only constrained 
by its coupling strength to the pair of initial quarks and its mass, which are the two free 
parameters in our analysis. Moreover, it is also equally probable that the vector diquarks 
have generation dependent couplings following Eq. 2.1. Therefore the couplings (X 2U , X u , X v ) 
involved in Eq. 2.1 are completely arbitrary and can in principle be large. Note that most of 
them are tightly constrained by flavor physics as they might mediate light meson or hadron 
decays [3, 27]. Therefore the constraints on the interaction of the vector diquark with the 
lighter quarks (first and second generation) are much more stronger, which means that vector 
diquark production at the LHC can have different allowed interaction strengths depending 
on the initial quarks participating in the production. To make our analysis more general we 
therefore choose to present our results normalized to the coupling strength. Where applicable, 
we would also assume that we work in the minimal flavor violating (MFV) scheme [45] for the 
couplings involving both the left- and right-chiral quarks with the vector diquark. It is worth 
noting that these colored states do not have direct coupling to a pair of gluons and thus the 
production cross section for diquark is limited by the flux of the initial partons in the proton 
at the LHC. However large QCD corrections can significantly alter the rates and modify the 
existing constraints on the mass and interaction strengths of such colored states. In this work 
we have chosen to ignore any electroweak corrections as interactions of the vector diquark to 
the electroweak gauge bosons might be model dependent. 


3 Production cross section at next-to-leading order 


We shall work in the “narrow-width” approximation where we can write the cross section as 
a product of the on-shell production and decay of the vector diquark (Vo) in a particular 
channel (XX) as 

a (pp XX) = 0- (pp -> V D ) x ISAzAL) (3.1) 

Thus a (pp —> Vo) gives the cross section for the production of the diquark resonance. The 
leading order or Born contribution to the on-shell vector diquark production comes from 
quark-quark initial states. The relevant Feynman diagram is shown in Fig.l. For the diquark 
of mass Mo, the parton-level cross section at the LO is given by 



- T 


where 

A 2 ttNo 
a ° ~ 2N% 


(3.2) 


(3.3) 
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Figure 1: Feynman diagram at Born level for the process qq —>• Ve¬ 


in the above, s is the partonic center of mass energy, Nq = 3 is the color factor of the quarks 
and r = M^/s. It is useful to rewrite the LO cross section in n = 4 — 2e dimensional form 
as this n-dimensional result will be used in the NLO calculation. Thus Eq. 3.3 can be put in 
the form 


< 5-0 


(n — 2)nN£)\ 2 (ii 2 )n 2e 

iiVg 


(3.4) 


Here A(/i 2 ) represents the running coupling parameter and p, defines the scale introduced to 
make the coupling dimensionless. From here onwards we shall drop the various indices from 
the coupling parameter introduced in the Lagrangian 2.1. The corresponding hadronic cross 
section at colliders can be obtained by convoluting the parton-level cross section with the 
parton distribution functions (PDF) of the initial quarks participating in the production, i.e. 


olo = — (q® q){ro), (3.5) 

where s is the hadronic center of mass energy and To = M^/s. We have used the notation 
for convolution of two functions, defined by 

(/i ® h)(xo) = [ dxi f dx 2 5(x 1 x 2 -xo)fi(xi)f 2 (x2). (3.6) 

io Jo 

Although the LO process involves colored particles only, the interaction strength does 
not involve the strong coupling g s but only the coupling strengths given by the free parameter 

o 2 

A. Therefore the one-loop QCD corrections at NLO are in leading order of a s = The 
0(a s ) QCD correction to the vector diquark production involves : 


• Virtual corrections due to one-loop gluon contributions. 

• Real corrections due to the gluon emission from initial state quarks and final state 
diquark. 

• For the complete 0(a s ) correction, one also needs to consider quark-gluon initiated 
diquark production with a jet. 
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We use dimensional regularization (DR) to regulate the ultraviolet (UV) and infrared (IR) 
singularities that may appear in these corrections. The renormalization of UV singularity 
and factorization of collinear singularity is carried out in the MS scheme. We have performed 
various checks, including the gauge invariance check with respect to the gluon at the amplitude 
and amplitude-squared levels, to ensure the correctness of our calculations. 



Figure 2: Feynman diagrams for virtual gluon correction to the process qq —>■ Vo- 


3.1 Virtual corrections 


The virtual corrections at 0(a s ) come from the interference of Born and one-loop amplitudes. 
The one-loop diagrams contributing to virtual corrections are displayed in Fig. 2. These 
diagrams are both UV and IR divergent. The required one-loop computation is carried out 
following the standard method of one-loop tensor reduction in n = 4—2e dimensions. We have 
listed all the one-loop scalar functions that we have used in the calculation, in Appendix B. 
The virtual cross section coming from vertex correction diagrams is given by 2 , 


OtsCe 


cry = cj b - 


2tt 


r r 8 i 

2vr 2 

77 i 

„ r i 2 

4 n 11 



+ TW 

+ Cf\ 2 

-+ 7T" - 8 \ 

L 13 euv 

3 

18 J 

1 £ uv e IR 

eiR > 1 


(3.7) 


The overall factor C e = r(i-e) (~uF~) a PP ears i n a H one-loop integrals regulated in DR. 
Cf a nd Co are the eigenvalues of the quadratic Casimir operator of SU(3)c acting on the 
fundamental representation and on the diquark representation respectively. For both the 
sextet and antitriplet diquark, Cf = 4/3 while Co is 4/3 for the antitriplet and 10/3 for 
the sextet diquark. The effect of external leg corrections can be incorporated in the wave 
function renormalization of the quark and diquark fields. Thus one can conveniently express 
the sum of Born and virtual cross section to 0(a s ) as [46], 


o'B+v = {Zl) 2 Z 2 d B + cry. 

2 We can also use this result to extract the vertex renormalization constant, 



(3.8) 
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The wave function renormalization constants and Z® for quark and vector diquark fields 

are, 


Zl 

Z? 




(3.9) 

(3.10) 


Note that these renormalization constants are calculated for on-shell quark and diquark fields, 
therefore, the IR singularity also appears. In DR, both are one as euv = £ir = e. However, 
the above form is suitable for extracting the full UV singularity in virtual corrections. The 
sum of Born and virtual cross section thus becomes, 


°B+v = 


1 + 


cy s C e 

2ir 

+ C D 


[C F 


- 7 , -h 7T — 8 


“IR 


Dr 


rll 

1 

1 

2vr 2 

771 


l 3 

euv 

eiR 

3 

+ isJ 

J 


(3.11) 


To get rid of the UV divergence in the above, renormalization of the coupling parameter 
A is necessary which is equivalent to adding an UV counter term of the following form to 
Eq. 3.11, 

(4vr Y 


-UV _ - a s 

a C.T. — 


11/1 /i 2 \ 

In— J 

I 3 Veuv [i-F 


'2irr ( l -e)”“ ° ' — ' “- 2 ' (3 ' 12 > 

where [ir is the renormalization scale. Hence the UV renormalized parton-level cross section 
to 0(a s ) for the production of diquark from qq initial state is given by 


o’B+v+C.T. = <7 B 


1 + 


<y. s C e 


{Cf(--|-- 

1 v e iR £ IR 
.2 


+ 7T - 


V eiR 3 s 3 18//j 


(3.13) 


Note that the procedure of renormalization has introduced a scale dependence in the cross 
section which would help in reducing the overall scale dependence due to the running of the 
coupling. After regulating the UV divergence, we are left with IR divergences, part of which 
will be canceled (due to Kinoshita-Lee-Nauenberg (KLN) theorem [46]) once we take into 
account the real gluon emission contribution. It is important to note that the singularity 
structure of virtual cross section is the same in the scalar [27] and vector diquark cases. Just 
like the singular terms proportional to Cf, we find that the singular term proportional to 
C D is also universal. 

Note that the results of this section can be utilized to predict the one-loop running of the 
quark-quark-diquark coupling A. The one-loop beta function due to 0(a s ) QCD correction 
is therefore given by 


/3(A) = M 2 


dlnA 

d/r 2 


47r 



(3.14) 
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Figure 3: Running of quark-diquark coupling with respect to the renormalization scale (hr) 
to ©(ab¬ 


solving this, the running of the renormalized coupling parameter X(hr) follows 3 


A (hr) = A (Q 2 ) 


! ^ s (h 2 r )11 ,(»r\ 

1_ ^ry CDln Mj’ 


(3.15) 


where Q is a reference scale which we will identify with Mjj (mass of the vector diquark) 
and choose X(M^) = 1. It is worth pointing out that in contrast to the scalar diquark 
case, the one-loop running of the coupling in the vector diquark case depends on the diquark 
representation and therefore will behave differently for the antitriplet and the sextet. This 
is highlighted in Fig. 3 where we show how the coupling A varies as a function of the 
renormalization scale hr■ Note that we have chosen A (Q 2 ) = 1 for Q = Mr> = 1 TeV as 
a reference point which is just for illustration purposes only. The scale dependence for the 
antitriplet vector diquark coupling is found to be at ~ 6% for the hr range considered while 
that for the sextet turns out to be significantly higher at ~ 16% for the same variation in hr- 
This is due to the dependence of the one-loop beta function onCo which takes different values 
for the two cases. Note that the running of the coupling will bring in a scale dependence for 
the LO cross section of the diquark too, similar to that observed for QCD cross sections due 
to the running of the strong coupling constant a s . 

3 We would like to point out that in the expression of running coupling for the scalar diquark case, given 
in Eq. 4.4 of Ref. [27], the factor of Cf should also be multiplied in the 0(a s ) term. 
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3.2 Real Corrections: qq channel 


Next, we compute the contribution from the gluon bremsstrahlung radiated from initial state 
as well as final state to 0(a s ). The process for the real gluon emission is, 

Qiipi) + Qj(P 2 ) ->• g{k) + V D (pi +P 2 - k). 

The Feynman diagrams which contribute to the NLO level gluon emission process for diquark 



Figure 4: Feynman diagram at leading order for the process qq —>• g Vp>. 


production is given in Fig 4. The full 0(a s ) spin and color averaged squared-amplitude for the 
three different diagrams shown in Fig. 4 can be expressed in terms of Mandelstam variables 
(s, t, u) in n = 4 — 2e space-time dimensions and is given by, 




12 — 


(n - 2)N D g 2 s X 2 f& 

4 Ngtu(t + u ) 2 

x ^4s 2 + (n — 2 )f 2 + 2 (n — 4 )tu + (n — 2 )u 2 + 4s(t + u)~j 


^ — Cntu + Cp(t + 


(3.16) 


where s = (p± + P 2 ) 2 ■, t = (pi — k ) 2 and u = (j >2 — k ) 2 . The partonic cross section for the 
real gluon emission process is obtained by performing the phase space integration in 4 — 2e 
dimensions and is given by 


a, 


R _ a 0 a s £ 


QQ 


s 2i r 




+ 4(1 + r 2 ) 


ir 

ln(l — t) 
1 — r 


+ C D \(— + ^)S(l-r)-l 

eiR 3 3 


1 + T 5 


CIR V 1 - T 


1 + r + T- 

1 — T 


+ O(cir) 


(3.17) 


In the above expression, the terms with (..,) + are the plus functions. The plus function 
distribution is defined in Appendix C. The IR divergence of real emission process originates 
from the phase space region where the emitted gluon is soft (/co —> 0) and/or it is collinear 
to the quarks. Since r = 1 corresponds to threshold production of the vector diquark, the 
1/eiR singular terms proportional to 5(1 — r) are due to the gluon becoming soft. On the 
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other hand, the 1 /ej R term arises when this soft gluon is also collinear to any of the two 
initial state quarks. The remaining singular terms in Eq. 3.17 are due to the gluon becoming 
collinear to quarks. Since the vector diquark is massive, the gluon emitted from it cannot 
be collinear thus explaining the absence of collinear singularity in Cd part of the expression. 
As mentioned above, the IR soft singularities cancel between real and virtual correction to 
qq —>• Vd- Adding the two cross sections given by Eqs. 3.11 and 3.17, we get 


^B+V+C.T.+R = (Tb+v+C.T. + d, 





1 + T‘ 


— T 


1 — T 

,3 s 3 18 / v 


2 fl + T + T 1 


T) 3 


1 — T 


(3.18) 


where we are left only with the collinear divergence terms as expected. The collinear di¬ 
vergences can be finally removed by redefining the quark PDF’s. In the MS factorization 
scheme, the universal counter term for collinear singularity is 


a. 


C.T. 

m 


dp a s (47r) £ 
s 2 -jt T(1 — e) 


/ 1 , n 2 

(-^ l n— T 

V eiR 


2 P qq (r) 


(3.19) 


where P qq (r) = Cp I pi 7 ). 1 is the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 

splitting function (probability of quark splitting into a quark and a gluon) and [ip defines 
the factorization scale. The total parton level cross section in qq channel is finally given by, 


Oqq = ds+V+C.T+R + d, 


C.T. 


s 


8(1-7-) + | 2P qq (r)ln( +C F [4(1 + r 2 ) 


2 J ln(l — r) 


+ C D 


2 f 1 + T + T 2 

3 V 1 -r 


1 — r 


2-k z 

+ (— - 8)5(1-r) 


+ 




2tH 


m d' 


+ 


143 


18 


)«5(1 


— T 


(3.20) 

The corresponding hadronic cross section is obtained by convoluting the parton level cross 
section with the initial state quark distribution functions, 


v qq — dr „ 

Jt 0 ' 


(q®q) 




qq- 


(3.21) 


If the initial state quarks are of different flavors q\ and qq, then replace, q®q —>• (qi®q 2 +Q 2 ®qi) 
in the above equation. 
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3.3 Real Corrections: qg channel 

As pointed out earlier, for a complete 0(a s ) contribution we should also consider the quark- 
gluon (qg) initiated process, 

Qi(Pi) + g(k) ->• V D (pi + k-p 2 ) + qj(p 2 ). 



Figure 5: Feynman diagrams for the process qg —» q Vo- 


The Feynman diagrams for this process are given in Fig 5. The total spin and color 
averaged amplitude-squared for the qg initiated process in terms of Mandelstam variables is 
given by 


V^TTTrTI _ 9 2 s >2p ie Np(Cpsu - C F (s + u ) 2 ) 

2-^ 19 2 Nq(Nq — l)sit(s + u ) 2 

x ^(n — 2)(s 2 + u 2 ) + 4 t(u + s) + 2(n — 4)srt + 4t 2 ^ 


(3.22) 


Note that the spin average for the initial state gluon introduces a term dependent on the 
space-time dimension n and also has a different color averaging factor compared to the qq 
initiated process for real corrections. However, as expected the above expression does match 
with that for the qq case without the spin and color averaging, under the interchange t -H- s 
and an overall sign. This is because of the crossing symmetry between qq and qg processes. 
The extra -ve sign in qg case results when one fermion is moved from initial state to the final 
state. 

The parton level cross section for the qg initiated process is 


R _ (T 0 a s 
q9 s 2n e 


-((1 - rf + r^ + 3 + 2r-3T* + ((1 _ T? + T>(1 _ t) 


2em 


+ 


Cp 
2 C F 


— 14-hr - 2 r 2 4- 2(1 4- r)lnr 


(3.23) 


where do is given in Eq. 3.3. As shown above, the cross section has IR collinear divergence 
which we remove by factorization in MS scheme. The required counter term is given by, 


C.T. = d 0 a, (4rr) e / 1 /x 
qg s 2 tt T(1 — e) v £ir. 


+ ln T2-) P 9 9 (' r ) 


(3.24) 
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with P q g(T) = i (1 - r) 2 + r 2 


. Hence the parton level cross section for the vector diquark 


production in qg initiated channel is given by 
T gg 


- _ ~R , C.T. 

- °99 + a qg 


<T() 

S 27T 


Pqg{ T 

A'jfr) 


1 \ g,p i / 


3 + 2t — 3r 2 


(3.25) 


+ 7 ^r"{ - 1 + ^ + t - 2t 2 + 2(1 + r)lnr| 


The corresponding hadronic cross section is obtained by convoluting the above parton level 
cross section with the initial state quark and gluon distribution functions, 


a qg ~ 


[ dr ^ [(<?< 

Jr 0 T* L 


, /To 

> g + g® fin¬ 


er, 


qg- 


(3.26) 


4 Decay Width: 0(a s ) correction 


Note that just like the LO cross sections for the production of the vector diquarks, the 
LO predictions for decay width of the particle also suffer from the renormalization scale 
uncertainties. Therefore for the sake of completion we would also like to estimate the effect 
of the QCD corrections on the decay width (T) of the vector diquark. Note that a primary 
requirement in assuming the narrow width approximation, one expects that the ratio T/Mp 
is relatively small and not exceeding ~ 10%. In order to remain in that regime, it is necessary 
to check that the decay width does not change by much under higher-order corrections. In 
this section, we compute the NLO QCD corrections to the vector diquark decaying into a 
pair of light jets, 

Vd(q) -1 Qi(Pi) + qj(P2)- (4.1) 

The leading order total decay width is given by 

F » = E < 4 - 2 > 


where i is the number of light quark generations which can couple to the vector diquark of 
a given electric charge. We have assumed that we can neglect all quark masses in the decay 
products (including top quark). The virtual corrections to the decay width involve the same 
Feynman graphs shown in Fig. 2 and has the same singular structure as given in Eq. 3.13 
for the on-shell vector diquark production. The same procedure followed in calculating the 
virtual corrections for the production cross section leads us to the UV renormalized virtual 
correction to the decay width which is given by 


T v 


= r n 


Us 

2n 


a 



+C D 


1 11 , 

-1-In 

{JR 3 


2 

3 
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Figure 6: Dependence of decay width on Mo (left) and on renormalization scale hr (right) 
to 0(a s ). To show the hr dependence we chose A {Mo) = 1 where Mo = 1 TeV. 


However we must point out that the real gluon correction is inherently different from that of 
the production. To compute the real gluon correction to the decay width, we need to consider 
the following three body final state, 


Vd( q) ->• QiiPi) + Qj(P 2 )+9- 


(4.4) 


Note that the calculation of real correction to diquark decay width requires three body phase 
space integration to be performed in n = 4 — 2e dimensions. For that we have followed the 
method given in Ref. [47]. The final expression with the real correction to the decay width 
is then given by, 


r =r„-c. 


„( 2 

3 

19 2^1 

( 1 

n\l 

C F ~ 2 ~ 

+ - + 

-X--7T 2 

+ Co -h 


\ e IR 

(IR 

2 ) 

\£IR 

3 J\ 


(4.5) 


By adding the virtual and real corrections to decay width all the singularities cancel as 
expected by KLN theorem. Thus, the complete NLO QCD correction to the diquark decay 
width is given by (from Eq. 4.3 and Eq. 4.5), 


Tnlo = Tq 
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(4.6) 


The corresponding expression for the case of scalar diquark is given in Appendix D. From Eq. 
4.6, we observe that the coefficient of C F is similar to that in SM (NLO QCD correction of an 
electroweak vector boson decaying into quark-antiquark pair) although here the final state is a 
quark-quark pair. We also find that a non-trivial contribution to the NLO decay width arises 
from the other Casimir, Co which takes different values for the two color representations of 
the vector diquark. 

We calculate the relevant A'-factor defined as the ratio of the NLO width to that of the 
LO width and plot it in Fig. 6. In the left panel of Fig. 6, we have shown the dependence of 
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the NLO K— factor for the decay width on the diquark mass. As /iR = Mr we can clearly 
see that the logarithmic term in Eq. 4.6 will not contribute and we should expect a constant 
value for a particular diquark representation. We however observe a slight variation for the 
NLO AT-factor for the widths of the antitriplet and sextet vector diquarks as we vary the 
mass Mr, which is only arising because of the running of the strong coupling a s (we have 
taken a s {Mz ) = 0.1184 as the reference value). We find that A'-factor for the sextet case is 
larger than the antitriplet due to larger Cr and increases the LO width by about 8 — 10% 
for the mass range Mr = 0.5 — 3 TeV. The corresponding LO width for the antitriplet vector 
diquark is modified less and increases by about 4.5 — 6% with the A'-factor. On the right 
panel of Fig. 6, we show the scale dependence ar of the decay width at LO and NLO and for 
sextet and antitriplet vector diquark states. As a reference point, we have chosen A(Afo) = 1 
where Mr = 1 TeV and we vary ar between Mr/ 2 to 3 Mr. The LO scale dependence is 
entirely due to the running of the coupling (see Fig. 3). We can clearly see that the inclusion 
of 0(a <J correction has significantly reduced the scale dependence. As one would expect, 
due to the smaller color factor Cr the scale variation for the antitriplet case is also smaller 
as compared to the scale variation for the sextet case. 

5 Numerical Analysis and Results 

In this section, we discuss the LO and NLO results for the vector diquark production at the 
LHC. We have used the CTEQ6L1 (CTEQ6M) [48] PDF’s for the parton fluxes in the colliding 
protons for our LO (NLO) results. In our calculations we choose /up = AR = Mr as the central 
scale for factorization and renormalization unless otherwise stated. Using our analytic results 
for the vector diquark production derived in the previous section we can now study how the 
cross sections are affected as a function of the collider center of mass energy (y/s) as well as 




Figure 7: Vector diquark production cross sections for the sextet and antitriplet cases at 
LO and NLO through the uu, uc and cc initial states as a function of the pp hadronic center 
of mass energy. 
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for different values of the mass (Mjj) of the vector diquark. The LHC has already completed 
its run at two different yfs of 7 and 8 TeV and there are plans of running the machine at 13 
and 14 TeV while future upgrades to 33 TeV is also possible. In Fig. 7 we show the LO and 
NLO hadronic cross sections for the on-shell vector diquark production as a function of the 
proton-proton collider center of mass energy, for a fixed value of Mp = 1 TeV. Note that the 
variation observed in the LO cross section can be attributed to the initial parton PDF’s only 
where, as the center of mass energy rises the on-shell condition of the diquark production for 
Me, = 1 TeV forces the colliding partons to carry a much smaller x (momentum fraction) 
of the proton beam energy. Therefore the initial quark’s flux grows giving rise to increase 
in the production cross section. The variation of the NLO cross section is however governed 
by both the partonic cross section and the PDF’s although the feature attributed to the LO 
behavior due to the PDF’s is similar. The plot is shown for three different quark-quark initial 
states, namely uu, cc and uc. It is worth recalling the fact that the coupling of the vector 




M D (GeV) 


Figure 8: Production cross section of the sextet and antitriplet vector diquarks at LO and 
NLO through the uu, uc and cc initial states as a function of the diquark mass Md at LHC 
with yfs = 13 TeV. 


diquark can be generation and flavor dependent. Therefore one can consider the diquark 
to be produced through initial partons of a particular fermion generation and flavor or it 
can be produced, mediated by interactions between different generations. We have chosen to 
normalize the cross sections with the coupling strength A squared so that it does not play 
a role here. Also note that although we always choose X(Mu) = 1 we have neglected the 
effect of the running of the coupling constant A in Fig. 7. Quite clearly, cross sections for the 
valence quark initiated processes are significantly large and reach appreciably high rates of 
above ~ 100 picobarns (pb) for 0(1) coupling strengths. Even the sea quark rates rise from 
a few 100 femtobarns (fb) to few ten’s of picobarns for both the sextet and antitriplet vector 
diquarks for 0(1) coupling strengths. When compared with the scalar diquark production 
rates we note that the LO cross section for the vector diquark production is exactly twice 
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that of the scalar diquark. 4 Again, as against the scalar case where same flavor initial states 
are disallowed for the antitriplet case because of the antisymmetric property of the K a ^, one 
gets all modes contributing in the vector case [35] . Thus a vector diquark which transforms as 
an antitriplet under SU(3)c would be produced through the initial valence uu and dd states 
resulting in a much higher cross section for the dijet final state compared to the scalar diquark 
which would have dominant production mode through ud initial states. One important point 
to note here is that if only flavor diagonal couplings are allowed for the uu type interactions 
then the vector antitriplet diquark will mediate same-sign top pair productions while the 
scalar diquarks will not, which would be a very interesting signal at the LHC. 
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Figure 9: Illustrating the NLO A-factors for the production of both sextet and antitriplet 
vector diquark at the LHC with y/s = 13 TeV, through the initial states uu, uc and cc as a 
function of the vector diquark mass M/> 

Since the vector diquark mass (Mq) is a free parameter, it is also instructive to know 
how the production cross section varies as a function of the diquark mass. We plot both 
the LO and NLO cross sections as a function of Mp at the LHC run with yfs = 13 TeV 
in Fig. 8. The plot is again shown for three different initial state combinations of quarks, 
namely uu, cc and uc. All these would lead to the production of a vector diquark of charge 
+4/3. The coupling strength has been factored out as before. We have varied Mjj in the 
range between 500 GeV to 1.5 TeV. Due to phase space suppression, the cross section goes 
down as we increase Md- It is worth pointing it out here that due to the difference in Njj, the 
sextet diquark production cross section at LO is just twice that of the antitriplet production 
cross section (see Eq. 3.3). However, the NLO cross sections are markedly different for the 
two cases and therefore the NLO cross sections for the sextet are no longer twice that of the 
antitriplet production. This will be evident from the A-factor estimates which we show next. 
Note that as all the different charged vector diquark productions are driven by the same color 

4 In Ref. [27] the interaction Lagrangian has an extra factor of 2^/2 thus giving overall rates higher than 
what we get here for the vector case. However once that is taken into consideration, one gets larger rates for 
the vector case as expected. 
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algebra for a given representation of SU(3)c the cross sections for them are eventually driven 
by the initial quark PDF’s that participate in the production. Therefore the nature of the 
plots for the production cross section for the |Q| = 2/3,l/3 charged diquarks is very similar. 



(t/M q (l/M 0 


Figure 10: Showing the scale dependence of LO and NLO production cross sections for 
sextet and antitriplet diquark states of mass Mjj = 1 TeV at the LHC with s/s = 13 TeV . 

In Fig. 9 we show the dependence of NLO A"-factor, defined as the ratio of the NLO 
cross section to the LO cross section, on the vector diquark mass Mo for both sextet and 
antitriplet diquark states. The A'-factors for the uu and dd initiated production are between 
1.5 and 1.3 for the mass range considered. We observe that the A'-factor for uu and uc initial 
states decrease with Mo while for cc initial state it increases which is mainly because of 
the difference in the PDF distributions for the valence and sea quarks in the proton. Also 
note that the AT-factors in the case of the vector sextet diquark are larger compared to their 
corresponding values in the vector triplet case which is unlike that observed for the scalar 
diquarks. For the scalar diquarks there is a partial cancellation between the Cp and Co 
terms, which gives a smaller AT-factor for the sextet case compared to the antitriplet [27], 
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Table 1: The LO and NLO cross sections (in pb) and K-factors for vector diquark production 
via different initial quark states at y/s = 8 TeV. We give the cross sections for both the sextet 
(S) and antitriplet (AT) diquarks. The uncertainties (in %) given for the cross sections are 
due to the the choice of scale Q = [i and is obtained by varying the scale from Mp>/2 to 2Mp>. 
We choose two reference values of the vector diquark mass Mp> = 1,3 TeV and a fixed value 
for the coupling, A = 1. 


while the Cp and Cp> terms in the vector case come with the same sign. However other 
features such as a larger K -factor for the sea quarks compared to the valence quarks remains 
the same, as this comes from their PDF behaviour as the factorization scale varies. 

One of the primary reasons for calculating the higher-order corrections to a scattering 
process is to minimize the scale dependence on measurable observables such as cross sections, 
that would affect the event rate estimates at experiments. We therefore make an estimate of 
the dependence of the choice of scale on the LO and NLO cross sections for the vector diquark 
production. To illustrate this we vary both the renormalization [ip and factorization [ip scale 
by a factor of two about the central scale /i = Mp> keeping [ip = [ip = /x throughout. Note 
that the renormalization scale dependence of the leading order cross section is governed by 
the one-loop running of the coupling parameter A. Thus the scale dependence of the LO cross 
section has an uncertainty of 0(a s ). Although, while predicting the scale dependence of NLO 
cross section, we should use two-loop running of the coupling, leading to an uncertainty of 
0(a 2 s )\ in absence of the two-loop result for running coupling we use Eq. 3.15 for predicting 
the renormalization scale dependence for both the LO and NLO cross sections for the vector 
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Table 2: The LO and NLO cross sections (in pb) and K-factors for vector diquark production 
via different initial quark states at y/s = 13 TeV. All other choices are similar to that in Table 
1. 


diquark production at LHC with y/s = 13 TeV. We plot our results in Fig. 10, where we 
can see clearly how the scale dependence of the NLO cross section is significantly reduced 
compared to the LO cross section. While the LO cross section varies between ~ ±30% 
for the vector sextet diquark for the three initial states uu, uc and cc as /x varies between 
Mu/2 to 2 Mu, the dependence is reduced to ~ ±10% for the NLO cross sections. For the 
antitriplet vector diquark, the dependence is relatively less compared to the sextet, of about 
~ ±(12 — 14)% for the LO cross sections which gets reduced to ~ ±4% for the NLO result. 
Notice that the scale uncertainty in antitriplet case is much smaller than that in sextet case 
which is reduced further when the NLO results are included. This is because of the Cu 
dependence (see Eq. 3.15 and Eq. 3.20) which is smaller for the antitriplet (Cu = 4/3) 
compared to the sextet (Cu = 10/3). 

Note that we have till now chosen to illustrate our results with figures for only the 4/3 
charged diquark production that couple to the first two generations of the fermions. But 
we should also note here that the vector diquarks with the 2/3 and 1/3 charge can have 
substantial rates only affected by the initial PDF’s of the contributing quarks. So to put 
the rate of production for the different vector diquarks in perspective we calculated all the 
modes that could contribute to its production and present the LO and NLO cross sections in 
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the relevant channels with scale uncertainties at y/s = 8 and y/s = 13 TeV run of LHC. To 
highlight the cross sections we have chosen two representative values of diquark mass Mjj = 1 
and 3 TeV and fixed the coupling A = 1. We show the cross sections for LHC with y/s = 8 
TeV and 13 TeV in Table 1 and Table 2 respectively. We assume that the couplings of vector 
diquarks mediating quarks of different generations is suppressed. So out of the 15 possible 
combinations we only consider 7 combinations with no inter-generation vertices. One can 
clearly see that the valence quark contributions dominate, with the uu and dd contributions 
being a few orders of magnitude higher than cc and ss respectively for Mjj = 1 TeV in Table 
1. For the 3 TeV diquark, the difference in orders is nearly doubled. A similar behavior 
is seen in Table 2. It is quite easy to understand that this happens due to the PDF’s of 
the quarks in consideration and the momentum fraction x of the initial proton that they 
carry. However the notable thing to consider is the fact that due to quite small production 
cross sections for the diquarks produced through second generation quarks, even with order 
1 coupling, the mass limits on them would be considerably weaker compared to the diquarks 
coupling to the first generation. As we have already determined a rough order of magnitude 
by which the cross sections differ for the first and second generation vector diquarks, it would 
give us a comparative idea of the limits on their coupling and mass from that derived for any 
one generation. We already have updated limits from dijet data by both ATLAS and CMS 
collaborations at the LHC [23, 26]. We use Ref. [26] of the CMS collaboration to derive 
the limits on the vector diquark mass and coupling. The CMS collaboration has given the 
upper bound on the cross sections for different resonant mass values which can be compared 
with the parton-level resonant production cross section (cr) times branching fraction (£>) in 



M D (TeV) M D (TeV) 

Figure 11: The constraints on the mass Mu and coupling A at 95% C.L. for the sextet and 
antitriplet vector diquark states at the LHC with y/s = 8 TeV using the LO and NLO cross 
sections. The values y/4n represents the perturbative limit for A while A = \J2Ait gives the 
upper bound on the coupling for T/Mu < 10%. 


- 20 - 













the narrow-width approximation using aBA, where A is an acceptance factor ru 0.6 [26]. 
We use this to derive limits for the vector diquark (both sextet and antitriplet) mass Md 
and its coupling A which interacts only with the first generation quarks. As these would 
be contributions coming through the valence quarks with the largest rates, the limits on the 
diquark coupling to the second and third generation quarks would be much weaker. In Fig. 11 
we show the 95% C.L. constraints on the mass and couplings of the vector diquark produced 
through uu, ud and dd fusion using the dijet data from Ref. [26] . The plots illustrate that 
all values of Md and A which are above the curves are ruled out by the CMS dijet data at 
95% C.L.. Note that we assume that the vector diquark couples to only one pair of quarks. 
We also show the perturbative limit of A = \/Ak in the plots, while A = \f2A-K gives the 
upper bound on the coupling for T/Md < 10%. As expected the strongest limits are for the 
4/3 charged diquark which couples to uu. The NLO corrections do modify the constraints 
to give slightly stronger limits compared to the LO results. For example, given a fixed value 
of the coupling A = 0.5 we find the dd initiated LO result for the antitriplet vector diquark 
gives a lower bound of Md — 3.03 TeV whereas the NLO corrections improve the limit by 
about 100 GeV to Md — 3.12 TeV. The corresponding limits for the sextet vector diquark 
at LO ( Md — 3.32 TeV) changes to Md — 3.42 TeV at NLO. The corrections in the other 
modes are also found to be between ~ 50-100 GeV. We have chosen not to show the effect 
of the associated scale uncertainties on the limits obtained. It should suffice to mention that 
the bounds using the LO cross sections would incorporate a much larger uncertainty band in 
the constraints compared to the NLO which is evident from the details given in Table 1 and 
Table 2. Also note that as the cross section for the second generation induced productions 
are at least 2 or more orders of magnitude smaller for similar couplings, the limits on the 
couplings would be relaxed by a factor of 10 or larger, allowing larger couplings for similar 
diquark mass. However one clearly finds a large parameter region still allowed for vector 
diquarks which should be explored at the upcoming run of LHC with y/s = 13 TeV. 

6 Summary 

In this work we have calculated the NLO QCD corrections to the vector diquark production 
at hadron colliders, namely the LHC. As colored particles are surely to be produced with 
large cross sections at hadron colliders, the discovery of any such state could be the first 
step towards discovering BSM physics at the LHC. Colored particles such as the vector 
diquark can mediate larger production rates for dijet and multijet events. We show how the 
NLO corrections to the vector diquark production affects the cross sections for the sextet 
and antitriplet representations. As the vector diquark couplings to the quark pair can be 
generation dependent, we find that valence quark processes have A-factors in the range of 
1.5 to 1.3 for a mass range of 0.5-1.5 TeV which decrease as we go higher in mass. The 
sea quark initiated production modes are found to have increasing values of the AT-factor 
as the diquark mass is increased. We also find that unlike the scalar diquarks, the sextet 
vector diquark has larger NLO corrections compared to the antitriplet. We also illustrate the 
scale uncertainties in the cross section for both the sextet and antitriplet vector diquarks and 
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find that the sextet vector diquark exhibits bigger scale uncertainty at LO compared to the 
antitriplet. The NLO corrected cross sections for both cases are found to show much lesser 
dependence on the scale variation. We also calculate the NLO corrections to the width of the 
vector diquark decaying to a pair of quarks. As a narrow-width approximation is considered 
large corrections to the width can affect predictions for relevant final states. We find that 
the A"-factor for decay width of the sextet diquark is around 1.08 — 1.1 while it is around 
1.05 for the antitriplet which is relatively smaller than that for the production cross section. 
However the scale uncertainties are relatively large for the decay width which get reduced by 
taking the NLO corrected widths. 

We have calculated cross sections for the vector diquark production at LHC with yfs = 8 
and 13 TeV arising from different generation quarks. We use the dijet data from the CMS 
collaboration for LHC with y/s = 8 TeV to put limits on the vector diquark mass and its 
coupling. We find that a large parameter region is still allowed for vector diquarks which 
should be explored at the upcoming run of LHC at y/s = 13 TeV. The current limits by 
the LHC experiments on the resonant particles include scalar diquarks but do not include 
vector diquarks. We have shown that using the same data one could also search for the vector 
diquarks and give search limits for such particles. 
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A Feynman rules 


The interaction Lagrangians given in Eqns 2.1 and 2.2 give the following Feynman rules (all 
momenta incoming) : 

• <TS(Pl)q'b(P2)Vf(P3) : - 6 qq ,KlP r .) 

where P T (P T ') can be Pl/r{Pr/l)- 

• V[ tl (p 1 )Vf 2 (p 2 )G A ^(p 3 ) ; -ig s T£[g^( Pl -p 2 r* + g^(p 2 -p 3 ) w 

+ 5 M3M1 (p3 -p\Y 2 } 


B One-loop scalars 


Here we list various tadpole (Ho), bubble (Bq) and triangle (Co) scalar integrals required 
in the calculation of virtual corrections in sections 3 and 4. For simplicity we take out 
the universal one-loop factor from these integrals which arise in DR and use the following 
notation, 
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We have labeled the UV and IR singularities of scalar integrals explicitly in our calcula¬ 
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The derivative of bubble function in Eq. B.7 
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used in the calculation of Z (f and ZP. 
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C Plus function 


For a function f(x), singular at x = 1. and a smooth function g(x), the plus function is 
defined by the following relation, 


[ dxf + (x)g{x) = [ f(x)[g{x)-g( 1)]. 

Jo Jo 


(C.l) 


Few plus function related identities which have been very useful in the calculation of real 
corrections are, 


f dxf + {x)g(x) = ( dxf(x)\g(x) - g{l)\-g{\) j dxf{x ) 
J a J a JO 

"ln(l — r)' 


(1 - r)( 1+2e ) (1 — t ) + 


- 2e 


1 — r 


" 2e S{ ' ~ T) 


f(x) 

(1 - t ) 


+ 


/(* 


( 1 - 


hi -M(l-r) 

T) J+ Jo 1-2 


(C.2) 

(C.3) 

(C.4) 


D 0(a s ) Correction to scalar diquark decay width 


The NLO QCD correction to the decay width for scalar diquark decaying into a pair of light 
jets is given by, 


Tnlo = C 




1 2 3 


where, the LO decay width Tq, is given by 



(D.l) 


A 2 

To = 77— Md- 
I 07 r 


(D.2) 


Note that the Cf part is exactly the same as one gets in the NLO QCD calculation of H —>• bb 
decay width [49]. We have used the following interaction Lagrangian for the scalar diquark 
(<Fj) case, 


^ = 7 \^Q C aK b Prq b + h.c.] + - M 2 $t $i . (D.3) 

i 1 + d qq') 

It should be noted that the coupling of the scalar diquark with two same flavor quarks is zero 
in antitriplet case. 
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E Useful relations 


Some of the relations among color factors that we have used to simplify various expressions 
in sections 3 and 4, are given below. For a more complete list one may refer to Ref. [27]. 


tabtA,bc = C F 5 ac (E.l) 

TfiT A , jk = C D 5 ik (E.2) 

Tr(iUiU) = N d (E.3) 

T \{K i t A t A K i ) = C f N d (E.4) 

T \{KH A K\t A ) T ) = ±- 2 C F N c (E.5) 

T a Tt(KH A K 1 ) = T$Tr(K\t A ) T K j ) = \c D N D (E.6) 

^zC f Nq = — 2C f N f + C F N F . (E.7) 


In the above t A b are the SU(3)c generators in fundamental representation while T, A are the 
generators in the diquark representation of SU(3)c- 

To calculate the real corrections to the 2-body decay of the diquark, the following relation 
has been used in simplifying the three body phase space integration in n = 4 — 2e dimensions. 

T(2z) = ^J-T(z)t(z+±). (E.8) 
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